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ABSTRACT 



Context. Water together with O2 are important gas phase ingredients to cool dense gas in order to form stars. On dust grains, H2O is 
an important constituent of the icy mantle in which a complex chemistry is taking place, as revealed by hot core observations. The 
formation of water can occur on dust grain surfaces, and can impact gas phase composition. 

Aims. The formation of molecules such as OH, H2O, HO2, H2O2, as well as their deuterated forms and O2 and O3 is studied in order to 
assess how the chemistry varies in different astrophysical environments, and how the gas phase is affected by grain surface chemistry. 
Methods. We use Monte Carlo simulations to follow the formation of molecules on bare grains as well as the fraction of molecules 
released into the gas phase. We consider a surface reaction network, based on gas phase reactions, as well as UV photo-dissociation 
of the chemical species. 

Results. We show that grain surface chemistry has a strong impact on gas phase chemistry, and that this chemistry is very different 
for different dust grain temperatures. Low temperatures favor hydrogenation, while higher temperatures favor oxygenation. Also, UV 
photons dissociate the molecules on the surface, that can reform subsequently. The formation-destruction cycle increases the amount 
of species released into the gas phase. We also determine the time scales to form ices in diffuse and dense clouds, and show that ices 
are formed only in shielded environments, as supported by observations. 

Key words, dust, extinction - ISM: abundances - ISM: molecules - stars: formation 



> 1. Introduction 



There is water everywhere in the cosmos. In cold and dense re- 
gions of the interstellar medium, dust gr ains become covered 
I by icy mantles mainly made of water ( Williams et al.l 119921: 

■ iGibb et aP l2004l) . In the star formation process, the ices that 
' were retai ned onto dust grains evaporate and drive a very rich 

■ chemistry JCaselli et al.lll993l: ICazauxet a"Lll2003h . At a latter 
stage, when planets are formed, the water still trapped in solid 
form is carried by asteroids and comets, can be brought to the 
different planets, giving birth, on our own planet, to oceans. 

Gas-phase chemical theory predicts that H2O and O2, along 
with O and CO, constitute the gas-phase reservoirs of ele- 

■ mental oxygen in molecular gas that is well shielded from 
ultraviolet radiation (e.g., Millar 1990). As a consequence, 
water togeth er with O 2 can be important coolants fo r dense 
gas (Goldsmit h &Lang er 1978; Hollenbach 1988; Neufel d et all 
ri995h . 

The formation and constitution of icy maiitles has been stud- 
ied by several authors (iTielens & Hagenll982llCuppen & HerbstI 
|2007|) . Observations revealed that icy mantles are sensitive to 
hard radiati on, and therefore co uld only exist in shielded en- 
vironments dWhittet et al.l l200 ll Papoular 2005]). Ices are sup- 
posed to be the result of gas phase species that freeze out onto 
dust grains, that were previously formed on the grain or in the 
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gas phase. Because the binding energy of water on bare carbona- 
ceous grains is lower than wat er on ice, the fi rst layer of ice is es- 
sential to form the icy mantle (Papoular 2005) and could have an 
important impact on the snow line position (Marseille & Cazaux 
submitted). Note that f or water sticking on bare silicate grains 
this could be different ( iGoumans et aljl2009l[Schilll973h 

Because of the inefficiency of gas phase routes, grain 
surface chemistry is required to e xplain the observed abun- 
dances of water in mo lecular clouds ( IP'Hendecourt et aPl 19851 : 
iHasegawa et al.lfl992h . H2O abundances can range from < 10"^ 
in the coldest clouds, that can be explained by freeze-out of 
water on dust, to > 10 '^ in warm gas and shocks, due to ice 
evaporation or spu ttering and exothermic gas p hase reactions 
(iBerginetalJI 19981 2000. lHollenbach et alJl2009i) . 

In this work we will investigate the formation of molecules 
such as OH, H2O, H2O2 and their deuterated forms in differ- 
ent astrophysical environments. In particular, we will concen- 
trate on grain surface chemistry, and the amounts of species 
that this process releases into the gas phase. The interac- 
tions between species and carbon surfaces (gr aphite or amor- 
phous carbon) has rec e ived a lot of atten tion (Ghio et al. 198_ 
'Pirr oiiello et aP 119971: Ijeloaica & SidisI [T999.: .Sha & Jackso'i, 
20021: [Bergeron et alj|2008l) . Theoretical and experimental stud- 
ies of the interactions between chemical species and these types 
of surfaces represent a considerable number of constraints for 
our model. Silicates, on the other hand, are more complex sur- 
faces that are more difficult to treat theoretically, resulting in 
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fewer constraints for the interaction s between speci e s and sur- 
face. However, recent studies from 'Goumans et alJ (12009^ re- 
vealed that the formation of water differs strongly as it is formed 
on olivine or on carbon grains. On olivine surfaces, the formation 
of water occurs via species strongly bound to the surface that can 
dissipate the excess energy in an efficient way via phonons. This 
allows most of the newly formed species to stay on the surface. 
Therefore, its seems that the formation of species on olivine sur- 
faces would not impact the gas phase significantly. In this study, 
because we want to assess how the chemistry on dust grains can 
impact interstellar gas, we exclusively follow the chemistry that 
occurs on carbon grains. To do so, we use Monte Carlo simu- 
lations to describe the formation of chemical species on such 
surfaces. Also, we consider several environments with different 
physical conditions, and establish how chemical species form on 
dust as well as their subsequent released in the gas phase. Finally, 
we discuss the formation of icy mantles in atomic and molecular 
clouds, and show that our model reproduces the fact that grains 
are covered by ices only in shielded environments. 

2. Grain surface chemistry 

2.1. Formation of water and its deuterated forms witli Monte 
Carlo simulations 



In a previous paper (ICazaux et al.ll2008l) . we discussed the for- 
mation of H2, HD and D2 on surfaces that are typical to the in- 
terstellar medium. We found that the formation of molecules de- 
pends on the binding energy of atoms with the surface and on 
the barrier that atoms from the gas phase have to cross to be- 
come strongly bound to the surface. Indeed, there are two types 
of interactions between the atoms and the surface: a weak one, 
called physisorption (van der Waals interaction), and a strong 
one, called chemisorption (covalent bound). Atoms from the gas 
phase can access easily the physisorbed sites and become ph- 
ysisorbed atoms. These weakly bound atoms can scout the sur- 
face at very low dust temperatures, meeting each other to form 
molecules. At higher dust temperature, the physisorbed atoms 
evaporate and the formation of molecules is insured by the con- 
tribution of chemisorbed atoms. The atoms from the gas phase 
can also enter in a chemisorbed site (and have to cross a barrier 
to do so), and can meet an already chemisorbed atom. This pro- 
cess will allow molecules to form for a wide range of dust grain 
temperatures (up to few hundreds of kelvins). 

We use a step by step Monte Carlo simulations to follow the 
chemistry occurring on dust grains. The dust grains are divided 
into square lattices of 100x100 adsorption sites. Each site on the 
grid corresponds to a physisorbed site and a chemisorbed site, 
so that the grain can be seen as 2 superimposed grids. Species 
that are coming from the gas phase arrive at a random time and 
location on the dust surface. This arrival time depends on the rate 
at which gas species collide with the grain. This rate of accretion 
can be written as: 



Race = nxVxo-S. 



(1) 



where n^ and Vx are the densities and velocities of the species 
X, cr is the cross section of the dust particule and S is the 
sticking coefficient of the species with the dust. We consider 
S=l, meaning that when a species arrives at a point of the grid, 
it can become physisorbed, or, if its chemisorption states ex- 
ists and its energy is high enough to cross the barrier against 
chemisorption, becomes chemisorbed. Jones & Williams (1985) 
show that the sticking coefficient depends on the proportion ph- 
ysisorbed/chemisorbed sites on the surface, and that for a gas 



at lOOK, most of the hydrogen atoms populate the chemisorbed 
sites. In our model, on the other hand, because of the high bar- 
rier to access chemisorbed sites, H atoms mostly arrive from the 
gas phase in physisorbed sites. In this sense, we overestimate the 
sticking coefficient since all atoms that aiTive on the grain stick. 
The species that are present on the surface can go back to the 
gas phase because they evaporate. This evaporation rate of the 
species i can be written as: 



Rev, 



'ap(i) 



Ei 
v/Xexp(-— -) 



(2) 



where E, is the binding energy of the considered species i, as 
reported in Table [1] 

The species that arrive at a location on the surface can move 
randomly, on the surface through tunneling effects and ther- 
mal hopping. The diffusion rates R,y for an atom to go from a 
site i (physisorbed or chemisorbed) to a site j (physisorbed or 
chemisorbed) can be written as: 



Rij = V; X Pij. 



(3) 



where v,- is the oscillation factor of the atom in the site i (which 
is of 10'^ s"' for physisorbed atoms), and P,j is the probabil- 
ity for the atom to go from a site i to j by tunneling effec t 
or thermal hopping, as described in ICazaux & TielensI (l2004l) . 
When thermal hopping dominates, which is usually the case for 
species going from a physisorbed sites to another physisorbed 
site, or for high temperatures, this probability can be written as 

E 

Pij - exp(-p^), where Eij is the energy of the baiTier between 
the site i andj, and T the temperature of the surface. The hopping 
baiTiers to go from a physisorbed site to another physisorbed site 
are taken as 2/3 of the binding energy (E,, listed in tablelTJ of the 
species on the surface, with a width of 3 A, which is what we de- 
rived for the case of H2 for mation of carbonaceous and silicate 
surfaces (ICazaux & Tielensll2004) . The baiTier for a physisorbed 
atom to become chemisorbed is discussed in section 2.2.2. 

Species present on the surface can also receive a UV photon 
and become dissociated with rates reported in section 2.4. Once 
2 species meet in the same site, they can form a new species if the 
activation barrier for the formation can be overcome. The dif- 
ferent reactions that can occur on the dust grains as well as their 
activation barriers £„ are listed in Table [T] When the barrier for 
a reaction has not been studied on the surface, we assume that 
the reacti ons are similar to gas phase reactions and use the NIST 
database dManion et al.ll2008h . However, recent studies showed 
that reactions that have a ba rrier in the gas pha se can be bar- 
rierless on the dust surface (lloppolo et al.ll200o ). Each time a 
reaction with a barrier occurs, we calculate the rate for crossing 
the barrier, which is Re^, = v x P^^, where v is the oscillation 
factor of the considered species, and Pe^, is the probability for a 
chemical species thermalized on the grain at a temperature Tdust 
to cross the activation baiTier £„. We assume here that these bar- 
riers have a width of lA. We compare the rate of crossing the 
baiTier Re„ to the rate for the chemical species to move out of 
the site R,y. If the rate of leaving the site is lower, then the re- 
action will occur, if not, the species will leave the site and go 
to a neighboring site. In this study we assume that species that 
remain on the grain are thermalized, and therefore that there is 
no preferential direction for the dissociated products that stay on 
the surface. We therefore consider that the probabilities to visit 
neighboring sites are equal. 

Once a chemical species present on the dust moves, or ad- 
sorbs a photon, or meets another species, the next event that will 
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occur to this species is determined as well as its time of occur- 
rence. Therefore, the events that concern every species on the 
dust are ordered by time of occurrence, and for each event that 
occurs, a next event for the concerned species is determined. 
The chemical species considered in our model are the follow- 
ing: H, D, O, OH, OD, O2, H2O, HDO, D2O, O3, HO2, H2O2, 
DO2, HDO2 and D2O2. Apart from H and D, that can also be- 
come chemisorbed as discussed in the next section, all the other 
species only fill the physisorbed sites. 

With the inclusion of deuterium, different products can be 
formed for the same reaction. We report the branching ratios of 
these reactions in table [T] For example, if we consider the reac- 
tion HD + O, the HD molecule has to be dissociated, and the 
oxygen reacts either with H or D. For this reaction, the branch- 
ing ratio are strongly dependent on the temperature. At low tem- 
peratures, tunneling dominates the reaction and HD + O gives 
mostly OH + D (with a rate kno) and very few OD + H Qs-dh)- 
At 300K the ratio ^ is of the order of 15 (iRobie et al.lli987l ), 

'<DH 

and strongly increases as temperature decreases. Therefore, we 
assume that at low temperatures, the branching ratio koH is neg- 
ligible compared to kHo- For the case of H2O2 + D, the reaction 
of one of the H atoms with a D atom, which gives HD + HO2, 
then breaks the 0=0 bound, and finally adds the free oxygen to 
HD. T he final products are therefore: OH + HDO (Albers et al. 
Il971h . For the other reactions, when no branching ratios have 
been reported (in the NIST database) we assume equal branch- 
ing ratios. 

2.2. Binding energies 

2.2.1. Physisorption 

The binding energies of the different species in physisorbed sites 
on car bon surfaces have been discussed by Cuppen & Herbst 
(l2007h and references therein. As discussed by these authors, 
the binding energies of species on the surface are poorly known. 
On amorphous carbon grains, it seems that only H and H2 
physisorption energies have been determined (Pii Tonello et al.l 
[1997, Katz et al. 1999). For most of the other binding energies, 
we refer to graphitic surfaces. H atoms are bourid to g raphitic 
surfaces with an energy of 550 K (^Berge ron et alJl2008l and ref- 
erences therein, Ohio et al. 1980), which is similar to the amor- 
phous carbon one. The physisorption energy of molecular hydro- 
gen on carbon nanoparticles (graphitic platelets) and polycyclic 
aromatic hydrocarbons (PAHs) is determined to be attractive be- 
tween 420 and 860 K, depending on the orientation of H2 and on 
the particle size. These energies are also d ependent on the tem- 
perature of the surface jHeine et al.l l2004). Therefore, we use a 
binding e nergy of 600 K, as i t has been determined fo r graphene 
surfaces (lAkai & Saitdl2003h . iBergeron et aP (l2008l) found that 
oxygen can be physisorbed on graphitic surfaces with energies 
of 1400 K, if the atom is placed on a top site (above a C atom) 
or on a bridge site (between two C atoms). This value is quite 
high c ompared to the valu e of 800 K used in ICup pen & Herbsn 
(l2007h , initially taken fr om lTielens & Hag en ( 1982). In a recent 
article, iLee et al] (l2009l) calculated the binding energy of ozone 
on graphene. Ozone molecules adsorb on graphene basal plane 
with binding energy of 2900 K, and the physisorbed molecule 
can chemically react with graphene to form an epoxide group 
and an oxygen molecule. The activation energy bartier from ph- 
ysisorption to chemisorption is very high (8000 K), and we will 
therefore neglect the fact that ozone molecules would overcome 
this bartier, and be dissocia ted. For the other s pecies , we con 
sider similar energies as in 



tea, tor tne otner s pecies , we c 
ICuppen & HerbstI (l2007h . The 



ergies of physisorption of the different species are reported in 
Table|2] 



2.2.2. Chemisorption 

Hydrogen atoms can chemisorb on graphite surfaces, but 
this process has an associated activation e n ergy of 0.2 eV 
("Ho rnekaer et"ai] l2006l ISha & Jacksoiil l2002l iJeloaica & SidisI 
199^. Recent studies show that once a first H atom is 
chemisorbed on a graphite surface, a second atom can chemisorb 
on the same benzene, in t he site opposite to this first atom (para 
site), without a bartier dHornekaer et al.l l2006l iRougeau et al.l 
1200 6*) and a third atom will form a molecule with one of these 
two adsorbed atoms without a barrier (Bachellerie et al. 2007'). 
This process is important to form H2 efficiently at high dust tem- 
peratures. This process will be added to our simulations in a 
coming paper (Cazaux & Morisset, in prep.). For now, we only 
consider that H atoms can become chemisorbed (with an energy 
of 0.7 eV) by crossing a bartier of 0.2 eV by tunneling effects 
(low temperatures) or thermal hopping (high temperatures), and 
approximate the baiTiers as being squared. Nevertheless, DFT 
calculations show that these baiTiers are far from being square 
resulting in higher rates for H atoms to become chemisorbed. 
This will be discussed in a future w ork (Cazaux & Morisse t 
in prep.) If the process discussed by iBachellerie et al.l (l2007h 
would be taken into account, once a first atom is present in a 
chemisorbed site, there is a barrier-less way to form and desorb 
a H2 molecule. In this study, we are thus underestimating the rate 
of H2 formation. This can have a small effect on the formation 
of molecules in diffuse clouds, which are mostly atomic envi- 
ronments, but has no effect for molecular environments, since 
the reactions to form molecules involve H2. 

Theoretical and experimental studies show that oxygen 
atoms chemisorb preferably above the middle of the C-C bond, 
creating a stable epoxide-like structure. Oxygen has adsorp- 
tion energies which are ~ 2.3-2.4 eV, with a bartier against 
chemisorption of the order of 0. 15-0.2 eV, depending on whether 
the O chemisorbs abo ve a C or in the middle of the C-C bond 
jFromherz et alj[l993h . In the singlet state, oxygen atoms can be 
chemisorbed with an energy of 2.5 eV (Jelea et al. 2004), and 
can react with H atoms to form OH radicals through an Eley- 
Rideal mechanism with an activation bartier that varies from 
0.17 eV to 0.29 eV. However, oxygen is mostly found in its 
triplet form in the ISM. Recent studies bv'Berge ronet alj (l2008h 
show that the chemisorption of oxygen in a triplet state has an 
associated baiTier of 0.2 eV, with a "meta-stable" chemisorbed 
state, meaning that the well lies above the energy level for ^O 
-I- graphite at infinite separation. Therefore, only for high tem- 
perature gas (2000 K) the chemisorbed state of oxygen (triplet) 
can be populated. In this work, because oxygen is mostly found 
under the triplet form in the ISM, we will consider that oxygen 
atoms cannot be chemisorbed. 



2.3. Desorption upon formation 

Once a reaction occurs on the surface, there is a certain prob- 
ability that the product of the reaction desorbs into the gas 
phase. This process depends on the energy released upon for- 
mation (enthalpy of reaction) and on the binding energy of the 
final product. The desorption upon formation concerns only ph- 
ysisorbed species that react with each other to form a product 
that is also physisorbed. Because of the excess energy contained 
in the molecule, when reactions have high enthalpies, this ph- 
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ysisorbed product can be ejected into the gas phase. This pro - 
cess has been discussed in a recent work bv lGarrod et alJ (l2007h . 
These authors quantify the probability of desorption using Rice- 
Ramsperger-Kessel (RRK) theory. However, this theory is de- 
signed for large molecules, and we expect a different behavior 
for diatomic molecules. With this method, they derive a prob- 
ability P that for a certain reaction, a certain energy E, higher 
than the binding energy of the product (E > Eg), is found in the 
bond. This probability can be written as: 



(4) 



£i-e 



where the energy released during the reaction is Zireac, and s is the 
number of atoms in the molecule which is produced. Once this 
probability is known, these authors use a method that depends 
on the energy loss on the surface, and that calculates the fraction 
of newly formed molecules that would desorb from the surface 
upon formation. This fraction is written as: 



vP 



aP 



v, + vP \ + aP 



(5) 



where Vj is the rate at which the total energy is lost to the surface, 
V is the oscillation factor of the molecule in the potential well, 
and a-—. We can speculate that the energy loss will depend on 
the porosity of the surface, meaning that a newly formed species 
can loose more efficiently its energy on a porous surface, since it 
will encounter more obstacles to bump into. Because water ices 
are supposed to be more porous than bare surfaces, molecules 
produced on bare grains should have higher probab ilities to des- 
orb di rectly into the gas phase upon formation iGarrod et aTl 
(l2007h use the data of iKroes & AnderssonI (l2005h to determine 
the rate Vj, and derive a factor a - 0.012. I n the case of bare 
grains , on the other han d, Katz et all (Il999l) . ICazaux & TielensI 



20041). Cuppen & Herbst (.20051 derive from the experiments of 
iPirronello et all (llQQTlT that the fraction of H2 that directly des- 
orbs into the gas phase upon formation is about 70 % for silicates 
and 60 % for amorphous carbon grains. With the theory devel- 
oped in Garrod et al. (2007), we obtain 1 %, which is far from 
what has been derived for H2. If we consider that the energy loss 
is much less efficient in the case of bare grains, meaning that the 
factor a has to be a few orders of magnitude higher than the 0.012 
determined by Garrod et al. (2007), we obtain a very similar frac- 
tion of desorption upon formation for H2 and H2O (if a = 0.7, we 
obtain 40 % and 39 % for the desorption upo n formation of H? 
and H 2O, respectively). In the discussion of ICuppen & Herbsn 
(1200 7"). it seems that the fraction of water t hat de sorbs upon for- 
mation is quite small (Kroes & Anderssonll2005h . Actually, this 
study deals with water ice, but for non-porous surfaces such as 

bare grains this fraction is differ ent. 

Recent studies from Berger on et al.l (l2008h . show that once 
physisorbed H and O atoms recombine on graphite surface, the 
nascent molecules can either directly desorb, if the translational 
energy is large enough, or the molecule is trapped in a quasi sta- 
tionary state where it oscillates in the physisorption well with 
a signicant amount of vibrational and rotational energy. As the 
H atom repeatedly collide with O, OH escap es from the sur- 
face (t he so-called complex mechanism). Also. Bachellerie et al.l 
(l2007h show that H2 molecules mostly desorb from the surface 
upon formation. This study shows that on a graphene surface at 
low temperatures, only 3% of the newly formed H2 molecules 
stay on the surface (they consider species staying more than 
12 ps as trapped). Therefore, we can assume that a significant 



fraction of newly formed species can desorb directly into the gas 
phase upon formation. 

We therefore will use an approximation that has been derived 
by Vasta et al. (2005, unpublished), where the desorption of H2 is 
high, and of H2O is low, depending on the binding energy of the 
final product and the enthalpy of reaction. This approximation, 
which reproduces the fraction of H2 that desorbs upon formation 
on silicates and amorphous carbon, also keeps this fraction small 
for water, and is more adequate for our purpose of modeling the 
formation of chemical species on bare grains. 



4.5x11600 



X (0.756 - 3 X IQ-^Es). 



(6) 



This empirical formula reproduces that 70 % (60 %) of H2 
molecules desorb upon formation on silicate (amorphous car- 
bon) surfaces, where H2 molecules have a binding energy of 
300 K (500 K). The fraction of molecules that desorbs upon 
formation linearly scales with the enthalpy of reaction. In our 
model, 60 % of the newly formed H2 molecules desorb upon 
formation, whereas this fraction is 36 % for the reaction H -i- O 
-^ OH and 15 % for OH + H ^ H2O. This fraction depends 
on the route to form the chemical species, since it becomes only 
4 % when OH is formed through O H- H2 and 0.8 % when H2O 
is formed through OH + H2. Therefore, the fraction of species 
released into the gas phase will vary strongly if the species are 
formed in atomic or molecular environments. 

2.4. Photo-dissociation and plioto-desorption 

Interstellar dust grains can be present in environments that are 
subject to radiation. Stars in the neighborhood can emit FUV 
photons that impinge on the interstellar dust grains. The flux of 
FUV photons can be written as follows (HoUenbach et al. 2008): 



Fpuv^GoFoexp-'-^^^ 



(7) 



with Fq ~10** photons cm"^ s ', the approximate local interstel- 
lar flux between 6 and 13.6 eV, Go a scaling factor for which the 
value 1 corresponds to the Milky Way interstellar radiation field 
in the FUV band, and Ay, the visual extinction. 

Once a UV photon arrives on a species, this species can 
be photo-dissociated. The products of the dissociation can be 
released directly into the gas ph ase, or one or both frag- 
ments can be trapped on the surface. lAndersson & van Dishoeckl 
(l2008h calculated the outcomes of the photodissociation of wa- 
ter molecules in water ices. Once a water molecule is broken 
into OH and H in the first monolayer of ice, several processes 
can occur with different probabilities if the ices are crystalline 
or amorphous ( Andersson et al. 2006) : a) H atoms are released 
into the gas phase and OH molecules are trapped (amorphous: 
92 %, crystalline: 70 %); b) H and OH are trapped (amorphous: 
5 %, crystalline: 14 %); c) OH and H desorb (amorphous: 2 % 
and crystalline: 6 %); d) OH and H recombine and form H2O 
which can be trapped (amorphous: 1 %, crystalline: 0.05 %) or 
desorb (amorphous:0.7 %, crystalline: 1 %). These results are 
different from crystalline and amorphous ices. In particular, H 
atoms due to dissociations are less trapped in crystalline ice than 
amorphous ice. According to these authors, the yield of desorp- 
tion of water molecules is rather small (0.2 %) for both typ es of 
ices, as confirmed by experiments from lOberg et al.l ( l2008h . 

In our simulations, we consider that an UV photon absorbed 
by a water molecule can lead to desorption of OH and H with 
a probability of 2 %. For the rest, we consider that the products 
stay on the grain surface as OH and H. Because we consider 
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bare surfaces, OH and H have different mobilities and there- 
fore the probability for the reaction and reformation of the wa- 
ter molecule should be different than on ices. We also consider 
that H atoms stay on the grain upon dissociation. Because of 
our lack of knowledge concerning the direct photo-desorption 
of species other than water, we only consider the desorption of 
water as OH and H into the gas phase. For the other species, 
we consider that there is no photo-desorption, but that photo- 
dissociation is possible, as discussed in the next paragraph. 
Th is seems to be a cor rect approximation, since, as discussed 
by lAndersson & van D ishoeck (2008), photo-dissociation (pro- 
cesses a and b discussed above) is more important than direct 
photo-desorption (process c). 



In this study, we consider that only the photons that arrive 
directly on the species can photo-dissociate them, meaning that 
the cross section of reaction is similar than the one in the gas 
phase. Table |3]presents the different photodissociation reactions 
and their associated photo-dissociation rates. Also, we consider 
that 2% of the products from a photo-dissociation reaction would 
directly desorb into the gas phase, while 98% of the same prod- 
ucts remain on the grain surface. We do not consider here the 
possibility of one product staying and the other one desorbing. 
We also take into account cosmic-ray-induced photons. These 
photons can interact with species present on the grain surface, 
and dissociate them with rates similar to the gas phase ones. 
These rates are also reported in Table [3] In our simulations, we 
consider that the products of the photodissociation are directly 
thermalized on the surface, meaning that they loose their excess 
energy in very short time scales. The 2 products of the dissocia- 
tion are placed in neighboring sites. 



2.5. Water clusters 



If a water molecule is located in a site next to another wa- 
ter molecule, a dimer can be crea ted, making the total binding 
energy larger. iBolina et all (l2005h and iBrown & Bohnal ilOOlh 
performed Temperature Programmed Desorption (TPD) exper- 
iments for the desorption of pure H2O ice on a graphite sur- 
face (highly ordered pyrolytic graphite). The experimental data 
showed that with increasing coverage, water molecules desorb 
at higher temperature. Water molecules are arranged in clusters, 
which have an increasing binding energy with increasing size. 
Therefore, as the coverage increases, the cluster becomes more 
important, and a higher surface temperature i s needed t o des- 
orb the molecules. Ab ini tio calculations from lLin et al.l (l2005h 
and lGonzalez et all (1200 7*). presented the adsorption energies of 
a single water molecule, but also a dimer, trimer, until the hex- 
amer. These binding energies increase with the number of water 
molecules that are present in the clu ster, but only if the cluster 
is in 2 dimensions. iLin et al.l (l2005h show that the binding en- 
ergy of a 3D cluster is mainly dependent on the number of water 
molecules which are close to the graphite surface. In this study, 
because we consider grain surfaces at low temperatures (<30 K), 
we do not take into account the formation of water clusters, since 
any newly formed water molecule that is not ejected upon for- 
mation stays on the surface. Water clusters and their impact on 
the freeze-out of water mantles and the snow line in proto-stellar 
environments is treated in another work (Marseille & Cazaux 
Submitted). 



3. Results 

3.1. Test case: the effect of the dust temperature and UV 
photons 

We calculate the number of molecules formed on the dust grain 
and released into the gas phase as a function of time for differ- 
ent dust grain temperatures. Note that we only consider grain 
surface chemistry, and that the gas phase reactions and abun- 
dances are not treated in this paper. We start our analysis with a 
high O/H and D/H ratio of 0.1, and «h = 10^ cm^^, in order to 
understand how grain temperature influences the chemistry oc- 
cuiTing on interstellar grains. We also add a UV radiation field 
to assess the changes in grain surface chemistry for the presence 
of photons. In this section we consider parameters that are not 
representative of any astrophysical environment in order to un- 
derstand how these parameters influence the chemistry on dust 
and its impact on gas phase. 

The population of the chemical species present on the dust 
grain and released into the gas phase, in absence of UV radia- 
tions, are presented in Fig. [T] When looking at the fraction of 
the species covering the grain (with monolayer equal to 1 cor- 
responding to a surface that is completely covered), it appears 
that the chemistry occurring on the dust grain is very different if 
the grain has a temperature of 10 K or 30 K. At 10 K, the atoms 
and molecules do not evaporate into the gas phase, and the H 
atoms (and D atoms) are the only species mobile on the surface. 
Therefore, the dominating process is hydrogenation and most of 
the oxygen is found in the form of H2O and HDO. On the grain 
surface, as shown in Fig.[Tl top left panel, H2O and HDO are the 
dominant species. Each time an OH species is formed and stays 
on the dust grain, a H (D) atom already present on the surface 
finds this OH (OD) radical and associates to form H2O (HDO). 
Therefore, the successive hydrogenations are very efficient, and 
O converts rapidly to OH and H2O. The species released into 
the gas phase, represented in fig.[T] top right panel, are identical 
to the ones formed on the grain: ~30 % of the oxygen coming 
on the grain is released in the gas phase as OH, and ~5 % as 
H2O. On 30 K grains, on the other hand, the H and D atoms 
evaporate, while oxygen atoms, more strongly physisorbed to 
the surface, stay and can encounter each other. Therefore, lots of 
oxygen bearing species are formed that contain more than one 
oxygen atom. The grain surface composition is rich in H2O and 
HDO (Fig. [1] bottom left panel). Lots of chemical species are 
only transitory, such as OH, OD, HO2 etc. and turn very fast into 
H2O, which explains their absence from the grain surface. The 
species released into the gas phase are more complex than grain 
surface species, as shown in Fig.[Tl bottom right panel. Oxygen 
atoms coming on the grain are released into the gas phase in the 
form of OH (20 %), O2 (10 %), H2O (5 %), OD (2 %) and H2O2 
(1 %). These species are released into the gas upon formation. At 
these temperature ranges, this process is the most important to 
populate the gas phase with species formed on dust since evapo- 
ration is negligible. 

When dust grains are subject to strong UV radiation fields, 
the species present on their surface can be photo-dissociated. If 
we consider 10 K and 30 K dust grains in an environment sub- 
ject to a UV radiation field of Go = 10^, and at an extinction 
of Ay - 1 mag, the molecules on the grain surface are photo- 
dissociated very fast (Fig.|2l left panels). Each time a molecule 
forms and stays on the surface, it can be destroyed by an inci- 
dent UV photon. Because of the high fraction of atoms on the 
surface, OH radicals are transitory species, and turn into H2O 
very fast. The fraction of H2O that covers the surface reaches an 
equilibrium between photo-dissociation into OH and H and ref- 
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ormation. The amount of water that is released in the gas phase 
upon formation is enhanced by the presence of the UV radiation 
field (Fig. |2] right panels). Without UV field, this amount is on 
the order of 5%. This fraction becomes much higher ~12-15 % 
in the presence of a strong radiation field. This is due to the fact 
that each water molecule that forms has a probability to be re- 
leased in the gas, but if this molecule stays, it can be dissociated 
and reform with the same probability of being desorbed. In this 
sense, the water molecules follow a cycle formation-destruction 
that enhanced the probability of being released in the gas phase. 
Analytically, we find that the amount of water that directly des- 
orbs into the gas upon formation should be ~5.4%(0.64x 0. 15 = 
0.054, amount of OH that stays on the surface upon formation 
times the amount of H2O that desorbs upon formation). Now 
when a strong UV radiation field is present, the molecules can 
be dissociated in products that can meet, and reform the same 
molecules. The probability that the H atom coming from the 
dissociation meet the other product OH to reform H2O is de- 
termined in the Appendix A. We determined that this process, 
commonly called backdiffusion, has a probability of 45% to oc- 
cur, and therefore enhances the fraction of water released into 
the gas phase by a factor 2 (~ 10% see Appendix A). This value 
is slightly lower than the one we found with our simulations, 
meaning that backdiffusion is enhanced by grain size limitation 
that we did not consider. Our calculations over the backdiffusion 
in App endix A differs from previous studies f rom Chang et al. 
( 120051) who studied the H -H H system and from Krud ( |2004|) who 
concentrated on stepped surfaces . In our study, only H atoms 
can move on the surface of bare grains, and we concentrated on 
the probability for H atoms resulting from the dissociation of 
water molecules to visit the location of the other product of the 
dissociation. 

It seems obvious that low dust temperatures favor hydro- 
genation and deuteration of the species on the surface. For higher 
dust temperatures, hydrogen and deuterium, can evaporate into 
the gas phase because of their low binding energies, while oxy- 
gen atoms, with higher binding energy, become mobile on the 
surface without evaporating. The resulting chemistry is oxy- 
genation which leads to a chemistry rich in oxygen bearing 
molecules. UV photons can dissociate most of the species on 
the surface. The products of the dissociation can meet again and 
reform the original species. Because each new species formed 
has a probability to be ejected into the gas phase upon forma- 
tion, the cycle formation-destruction-formation increases the to- 
tal fraction of species released into the gas phase. 

4. From diffuse to dense clouds and PDRs 

In this section, we discuss the formation of molecules on grains 
in different enviro nments. The parameters used in our simula- 
tions are based on ISnow & McCalll ( l2006h and are reported in 
Table |4] For each environments, we performed at least 4 differ- 
ent simulations with different seeds, and we report the errors bars 
that represent 95% of confidence level. 

We first compute the chemistry occurring on dust grains 
in diffuse clouds, where the gas phase composition is mostly 
atomic. We assume a density of 100 atoms cm""*, a visual extinc- 
tion of 0.5 mag and an impinging UV radiation field of G()=l. 
The species present on dust grains and released into the gas 
phase are reported in Fig. [3] In this environment, the process 
to form water is through the association of H and O atoms. The 
water formed on the grains can desorb or remain. The fraction 
of the incoming oxygen atoms that come and leave the grain as 
OH is ~ 35 %, while it is ~ 10 % as H2O. Indeed, an oxygen 



atom that arrives on a grain and meet a H atom, will form OH. 
Because of the energy released during this reaction, OH desorbs 
in 35 % of the cases. The fraction that remains on the surface 
(65 %) can meet another H atom to form H2O. Again, this re- 
action is exothermic and in 15 % of the cases, the product H2O 
is ejected into the gas phase, making a total of 10 % of the O 
transforming into gas phase H2O (0.65x0.15 = 0.1). The water 
that remains on the surface builds up very slowly until it reaches 
0.2% after ~ 100 years. This shows that almost no water ice cov- 
ers dust grains in D iffuse clouds, as supported by observations 
(IWhittetetalJl200ih . 

In translucent clouds, on the other hand, the medium is 
mostly molecular We assume a density of 500 atoms cm"-', a 
visual extinction of 1 mag and an impinging UV radiation field 
of Go = 1 . In this cloud, the deuterium and oxygen species are 
still in atomic form. Because in these clouds the hydrogen is in 
molecular form, the most efficient route to form H2O is H2 + O 
— >OH + H. At Tdust = 14 K, the hydrogen atoms evaporate, and 
water forms through the association of OH and H2. Because the 
enthalpies of reactions are much lower for the formation of H2O 
when the reaction involves H2 than when it involves H atoms, the 
desorption upon formation becomes 3.2% for H2 +0 and 0.88% 
for H2-1-OH. Therefore, the newly formed species mostly stay on 
the grain surface, and we find that no species are released into 
the gas phase ( see Fig.|4|. The first layer of ice is growing until 
it reaches a coverage of few percent, confirming the absence of 
ice on the surface of dust grains in translucent clouds. 

In dense clouds, the medium is also molecular, but denser. 
We consider a cloud with a density of 5000 cm^^, at a visual ex- 
tinction of 5 mag, and an impinging UV field of Go = 1 . At this 
extinction oxygen and deuterium are still in atomic form and the 
dust grains are colder, Tdust = 12 K. The formation of H2O in- 
volves the reaction H2 -H O — > OH + H. As discussed before, this 
reaction releases less energy than the reaction involving atomic 
H, and therefore only a small percentage of OH (~3 %) is re- 
leased into the gas phase. Because the dust grains are at lower 
temperature than in the translucent cloud case, H atoms do not 
evaporate. The products of this reaction, OH and H, that stay 
on the surface can react and form H2O. If these two products 
meet, 15 % of the formed H2O is released into the gas phase. 
Our results, presented in Fig.|5] show that more than 10 % of the 
oxygen is released as H2O into the gas phase. This means that 
the products H and OH have a probability of 60 % to meet each 
other The OH that does not meet the H atoms will encounter an 
H2 molecule and form H2O -i- H. The part released into the gas 
phase through this process is negligible. 

We also calculate the formation of molecules on dust grains 
in photo-dissociatio n regio ns. For this we use the results of 
'Me iierink & SpaansI ( l2005h for a PDR with a density of 10^ 
atoms cm~^, with a UV radiation field of Go=10^ and a visual 
extinction of 5 mag. Hydrogen is converted into H2 but oxy- 
gen and deuterium are still in atomic form. The grain and gas 
temperatures are as high as 30 K. For these grain temperatures, 
the chemistry is very different than from above, as shown in 
Fig. refenv4. H atoms as well as H2 molecules, because of their 
low binding energy, evaporate efficiently, while oxygen atoms 
remain on the grain to form O2 and O3 molecules. The forma- 
tion of water occurs through the successive hydrogenation of O2 
(Mi vauchi et alll2008l IIoddoIo et al.ll2008h . H +O2 ^ HO2 and 
HO2 -H H -> H2O2 and H2O2 h- H ^ H2O. An other important 
route is the hydrogenation of O3, as shown by iMokrane et al.l 
(,2009i) : O3-HH ^ OH -H O2 and OH -H H ^ H2O. 
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5. Discussion and conclusions 

We used Monte Carlo methods in order to follow the chem- 
istry involving H, D and O species on interstellar dust grain sur- 
faces. The formation of surface species as well as species re- 
leased into the gas phase is computed and allows us to estimate 
the efficiency for the formation of O bearing species. Our model 
shows a strong chemical differentiation with grain temperature. 
At Tdust^lO K, processes involving hydrogenation are favored, 
while rdust=30 K favors oxygenation. UV photons can dissociate 
the species present on the surface, and the product of this dissoci- 
ation can meet again and reform the initial species. Because lots 
of chemical reactions used in this work are exothermic, a fraction 
of the newly formed species desorbs directly into the gas phase. 
Therefore, the cycle formation-dissociation-formation increases 
the fraction of the species released into the gas phase signifi- 
cantly. 

We applied our model to different astrophysical environ- 
ments. In diffuse clouds, water is formed through the succes- 
sive hydrogenation of oxygen. These reactions are exothermic 
and allow an important fraction of newly formed OH and H2O 
to desorb into the gas phase (~35 % and ~15% respectively). 
In translucent clouds, on the other hand, the medium is molecu- 
lar and the formation of water involves molecular hydrogen and 
atomic oxygen. These reactions produce a very small amount of 
energy, and therefore no species are released into the gas phase 
(~ 3% for OH and ~ 1% for H2O). In diffuse clouds, the tem- 
perature of the dust grains is slightly lower, and H atoms do not 
evaporate. Therefore, as in the case of translucents cloud, oxy- 
gen meets H2 to give OH and H, but in this case the H atoms 
produced do stay on the grain. Water can form through the as- 
sociation of OH and this H (with a probability of 45%), and 
again, an important fraction of H2O is released into the gas phase 
(~7%). In environments with warmer dust grains, such as PDRs, 
the chemistry occuiTing on the grain and the products released 
into the gas phase are very diffe rent. The formation of water in - 
volves O2 (lloppoloetalj|2008l) and O3 dMokraneet alj|2009l) . 
and many oxygenated species are released into the gas phase. 

The formation of ices requires very long time scales 
in diffuse and transluc ent clouds, as akeady discussed by 
ICuppen & HerbstI (l2007h . In their work, these authors showed 
that a time of 10^ years is required to form the first monolayer 
of ice in dense clouds, while water coverage would reach a max- 
imum of 50 % of the surface in translucent clouds and 10 % of 
the surface in diffuse clouds. In our model, the water ice cover- 
age reaches 0.3 % in diffuse clouds and 1 % in translucent clouds 
after 30 years. We obtain a water ice coverage much lower than 
ICuppen & HerbstI ( 2007) since we consider that for each reaction 
that occurs on the surface, a certain fraction directly desorbs into 
the gas phase upon formation. In diffuse clouds, this fraction is 
of 36% and 15% for the formation of OH and H2O respectively. 
Also, our model differs from Cuppen & Herbst (2007) since it 
does not allow species to accrete on top of each other and there- 
fore does not consider the formation of ices. 

Even if our model does not allow us to estimate the coverage 
of ices and their morphology, the water coverage can be deter- 
mined in steady state as an equilibrium between water formation 
and destruction by UV radiation. The formation rate of water 
on dust is set by Rh,o= no vq ere, where no is the density of 
atomic oxygen in the gas phase, vq the thermal velocity of oxy- 
gen atoms (~4.6 10"* -J^cm s"') , cr the cross section of the 
grain and e the formation efficiency of water, that stays on the 
grain. This efficiency directly depends on the formation route 



of water in different environments. In the case of diffuse clouds, 
water forms through the reactions O H- H ^ OH, for which 64% 
stays on the dust grain (see section 2.3), and OH -1- H — > H2O 
for which 85% stays on the dust, yielding an efficiency on the 
order of 50 %. In translucent clouds, water forms through O + 
H2 ^ OH H- H (96% stays on dust) and OH H- H2 ^ H2O + 
H (99.2% stays on the dust), which brings the efficiency of wa- 
ter formation to ~95%. In dense clouds, this efficiency depends 
on the reactions O H- H2 ^ OH + H and OH + H-» H2O, and 
therefore is on the order of ~81%. Note that in the PDR case, 
water formation follows another route involving O2 that brings 
the efficiency to 5%. The destruction by UV can be written as: 
Rdes=8 lO"'" Go UhjO exp(-2.2Av). In steady state, the fraction 
of the water molecules that covers the surface of a dust grain can 
be written as: 



nHjO 



noVoO"sitee 



n^ite 8 X 10-"'Go exp(-2.2Av) 



= 1.72510" 



5 nHeyTg^ 
Go exp(-2.2Ay) 



(8) 



where cr^^xi^ is the cross section of one site (~10"'^cm^) and Uh 
the total density of hydrogen. This fraction is ~ 0.3% in dif- 
fuse clouds, which is consistent with our calculations that show 
that this coverage is reached after 100 years (time to reach the 
steady state; see Fig. |3]l. In translucents clouds, the density is 
higher, and the coverage at steady state is of the order of 5 %. 
For dense clouds, the medium is shielded, with an extinction of 5 
magnitudes. The coverage becomes larger than 1, meaning that 
dust grains are covered by icy mantles with several layers. In 
PDR environments, with an high impinging UV radiation field 
(lO^*), and for an extinction of 5 magnitudes, the coverage of ice 
is of 0.25 %. These estimates, confirmed by our Monte Carlo 
simulations, show that ices that cover dust grains can only ex- 
ist in shielded environments, if the UV radiation field is not too 
important. This is suppo rted by observation s of ices in diffuse 
and de nse environments dWhittet et alJl200l1) . lJones & Williamsl 
(|1984|) . estimate that the efficiency to form water on dust grains 
through repetitive hydrogenation of oxygen should be higher 
than 25% in order to reproduce the observations of the ice band 
in the Taurus molecular cloud. In this work, we find that the effi- 
ciency for the formation of water is at least of 50%, which agrees 
with their findings. However, we have established that the suc- 
cessive hydrogenation of oxygen to form water was not impor- 
tant in molecular clouds, but that the reaction of oxygen with H2 
is the dominant route. Indeed, even if an important barrier ex- 
ists for oxygen to react with H2, these reactions are possible on 
surfaces since the species can repetitively meet each other which 
increases the cumulative probability to overcome the barrier. In 
this sense, dust grain surfaces can provide natural place for reac- 
tions with high baiTiers to occur, while these reactions would be 
highly unlikely in the gas phase. 

A major result of this paper is that chemistry on bare grains 
has an important impact on gas phase chemistry because of 
the exo-thermicity of certain reactions on the surface. Once the 
grains become covered by icy mantles, the binding energies are 
much higher and the newly formed species remain trapped on 
the surface. Therefore, it would be very interesting to follow 
the chemistry that occurs in a diffuse cloud as it evolves into 
a dense cloud under the influence of hydrodynamics and gravity, 
and measure the impact of the grain surface chemistry on the gas 
phase composition while the grains are still bare grains. 

We should mention that the probability for chemical species 
to desorb from grain surfaces upon formation has been seen 
experimentally (Pirronello et al. 1993) and has been show n 
theoretically (.Bachellerie et aL ,20071 iBergeron et"an l2008h . 
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However, the dependance of this probability on the exo- 
thermicity of the reaction is unknown and has only been esti- 
mated empirically. The impact of grain surface chemistry on gas 
phase has been estimated in this work with the assumption that 
species forming on dust surfaces have higher probabilities to be 
released in the gas as the reaction is exothermic. Further studies 
are highly needed to assess this point and better understand the 
effect of grain surface chemistry on the gas phase composition. 
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Table 1. Reactions adopted in our chemical model 



Reaction 




Branching ratio 


Barrier in K 


references" 


H+H 


^H2 









H+D 


^HD 









H+O 


^OH 









H+OH 


^H20 









H+OD 


^HDO 









H+02 


^H02 




250 


Walchetal.('1988") 


H+H20 


-^ H2+OH 




9600 




H+HDO 


^ H2+0D 


50% 


9600 






^ HD+OH 


50% 


9600 




H+D20 


-^ HD+OD 




9600 




H+03 


^ 02 + OH 




450 


Lee etal. (•1978) 


H+H02 


^H202 









H+D02 


^HD02 









H+H202 


^ OH + H20 




~0 


foppolo et al. (2008) 


H+HD02 


-> OH + HDO 


50% 


~0 


loppolo et al. (2008J 




^ OD + H20 
-^ OH + D20 


50% 
50% 


~0 
~0 


IoDDoloetal.r2008) 


H+D202 


loDDolo et al. (2008) 




^ OD + HDO 


50% 


~0 


loDDolo et al. (2008) 


D+D 


^D2 









D+O 


^OD 









D+OH 


^HDO 









D+OD 


^D20 









D+02 


^D02 




250 


Walch et al. (1988) 


D+H20 


-^ HD+OH 


50% 


9600 






^H2+OD 


50% 


9600 




D+HDO 


^ HD+OD 


50% 


9600 






-> D2+0H 


50% 


9600 




D+D20 


^ D2+0D 




9600 




D+03 


^ 02 + OD 




450 


Lee etal. (1978) 


D+H02 


^HD02 









D+D02 


-4 OD + OD 









D+H202 


^ OH + HDO 


100% 


~0 


Jpppolo et al. (2008) 




^ OD + H20 


0% 


~0 


foppolo et al. (2008) 


D+HD02 


^ OD + HDO 
^ OH + D20 


50% 
50% 


~0 
~0 


loppolo etal. (2008) 




loppolo et al. (2008) 


D+D202 


^ OD + D20 
^OH+H 




~0 
3000 


loppolo et al. (2008) 


0+H2 




O+HD 


^OD+H 


0%'' 


3000 






^OH+D 


100%'' 


3000 




0+D2 


^OD+D 




3000 




O+O 


^02 









O+OH 


-^02 + H 









O+OD 


^02 + D 









0+02 


^03 









0+H20 


^ OH + OH 




8500 




O+HDO 


^ OH + OD 




8500 




0+D20 


^ OD + OD 




8500 




0+03 


^02 + 02 




2000 




0+H02 


^ OH + 02 









0+D02 


-> OD + 02 









0+H202 


^ OH + H02 




2000 




0+HD02 


^ OD + H02 


50% 


2000 






-^ OH + D02 


50% 


2000 




0+D202 


^ OD +D 02 




2000 




H2+OH 


^ H20+H 




2600 


Schifr(197J) 


H2+OD 


^ HDO+H 


100%"^ 


2600 


'Schiff(1973) 




-^ H20+D 


0%" 


2600 


ScMff(1973) 


H2+02 


-^ OH+OH 




35000 




H2+ H02 


^ H202+H 




10000 




H2+ D02 


^ HD02+H 


100% 


10000 






^ H202+D 


0% 


10000 
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HD+OH 


^ HDO+H 


50% 


2600 


ScMffQQTJ) 




^ H20+D 


50% 


2600 


S^chiff(197J) 


HD+OD 


-^ HDO+D 


50% 


2600 


^hifr(1973) 




-^ D20+H 


50% 


2600 


Schiffri973^ 


HD+02 


^ OH+OD 




35000 




HD+ H02 


^ HD02 + H 


50% 


10000 






^ H202 + D 


50% 


10000 




HD+ D02 


-^ D202+H 


50% 


10000 






-^ HD02+D 


50% 


10000 




D2+OH 


-^ HDO+D 


100%"^ 


2600 


Schiff(1973) 




^ D20+H 


0%" 


2600 


Schifr(197J) 


D2+OD 


^ D20 + D 




2600 


Schiffri973) 


D2+02 


^ OD + OD 




35000 




D2+ H02 


-^ D202+H 


0% 


10000 






-^ HD02+D 


100% 


10000 




D2+ D02 


^ D202+D 




10000 




OH+OH 


^H202 









OH+OD 


^HD02 









0H+02 


^ H02+0 




26000 




0H+03 


-^ H02 + 02 




1000 




OH+ H02 


-^ H20 + 02 









OH+ D02 


^ HDO + 02 









OH+ H202 


^ H20 + H02 









OH+ HD02 


^ HDO + H02 


50% 









^ H20 + D02 


50% 







OH+ D202 


-^ HDO + D02 


100% 









^ D20 + H02 


0% 







0D+02 


-4 D02+0 




26000 




0D+03 


^ D02 + 02 




1000 




OD+ H02 


^ HDO + 02 









0D+ D02 


^ D20 + 02 









OD+ H202 


-^ HDO + H02 


100% 









-^ H20 + D02 


0% 







OD+ HD02 


-^ HDO + D02 


50% 









^ D20 + H02 


50% 







OD+ D202 


^ D20 + D02 









02 + H202 


^ H02 + H02 




20000 




02 + HD02 


-^ H02 + D02 




20000 




02 + D202 


-^ D02 + D02 




20000 




° INIST Chemical Kinetics Data! 


)ase database: Mat 


ion et al. ("2005 


) if no ref. is in 



S. Cazaux, V. Cobut, M. Marseille, M. Spaans, P. Caselli: water on bare grains 



11 



Table 2. Physisorption and chemisorption energies of the considered species. 
Species Physisorption Chemisorption 



H 


550 


8500 


D 


550 


8500 





1390 




H2 


600 




HD 


600 




D2 


600 




OH 


1360 




OD 


1360 




02 


1440 




H20 


2000 




HDO 


2000 




D20 


2000 




03 


2900 




H02 


2160 




D02 


2160 




H202 


2240 




HD02 


2240 




D202 


2240 





Table 3. Photodissociation reactions 



Reaction 


"photis'^) 


r 


CRo'' 


H2 + photon ^ H + H 


5.0(-ll)'' 


a 


0.05 


HD + photon ^ H + D 


5.0(-ll)'' 


a 


0.05 


D2 + photon ^ D + D 


5.0(-ll)'' 


1.8 


0.05 


OH + photon ^ O + H 


3.9(-10)^ 


2.24 


254.5 


OD + photon ^ O + D 


3.9(-10r 


2.24 


254.5 


H20 + photons OH + H 


8.0(-10r 


2.2 


485.5 


HDO + photon-* OD + H 


8.0(-10)" 


2.2 


485.5 


D20 + photons OD + D 


8.0(-10)^ 


2.2 


485.5 


02 + photon-^ O + O 


7.9(-10r 


2.13 


375.5 


03 + photons 02 + 


1.9(-9r 


1.85 


375.5 


H02 + photon ^ OH + O 


6.7(-10)" 


2.12 


375 


D02 + photon ^ OD + O 


6.7(-10r 


2.12 


375 


H202 + photon ^ OH + O 


9.5(-10r 


2.28 


750 


HD02 + photon ^ OD + O 


9.5(-10r 


2.28 


750 


D202 + photon ^ OH + O 


9.5(-10)" 


2.28 


750 


ffphot = Cpho, X exp(-yAv) (s" 


');CRpho,= 10 


-"• CRo 


(s->) 



' depends on the H2 and HD self-shielding; * van Dishoelc (1998); ' van Dishoelc (2006); '' Woodall et al. (2007) 



Table 4. Parameters used in our 


simulations 
















Environment 


hh 


hhi 


nuj 


hdi 


noi 


Adust 


T«as 


Go 


Av 


Diifuse Clouds 


100 cm-^ 


100 cm-^ 


cm-' 


0.002 cm-^ 


0.03 cm-' 


18 K 


100 K 


1 


0.5 mag. 


translucents clouds 


500 cm-3 


6cm-3 


247 cm-' 


0.01 cm-3 


0.15 cm-' 


14 K 


50 K 


1 


Imag. 


Dense Clouds 


5000 cm-3 


50 cm-3 


2742 cm-3 


0.1 cm-3 


1.5cm-3 


12 K 


20 K 


1 


5 mag. 


PDR 


1000 cm-3 


10 cm-3 


495 cm-3 


0.02 cm-3 


0.3 cm-3 


30 K 


30 K 


10' 


5 mag 
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grain. lOK 



0.01 - 



0.001 




time (yrs) 



0.1 



0.01 - 



0.001 



grain 30K 



lime (yrs) 



H,0 




HDO 



S.0.1 



0.01 




1 2 

time (yrs) 



Fig. 1. Chemical species present on a grain surface of 10 Kin monolayer (1 monolayer is a total coverage of the surface) as function 
of time (left panels). Fraction of newly formed oxygen bearing molecules f(0) that are released into the gas phase under the form 
of other species (right panels). Top: temperature of the grain surface is set to 10 K. Bottom: temperature of the grain surface is set 
to 30 K. 
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U .^J I 
















gi'ain 


lOK; 


Go 


10^' 










J-lgU 


_^^ 


/ 


J 




HUU 




J — '^ 




L 


^ 




>> 

s 




HD 




^ 


y 






OH 


0.001 




' r 


n 




r 




— 


- 


f , 






D,0 




\ 













2 4 

time (yrs) 



1 - 



o 



'0.1 



0.01 



gram lOK; Go=103 



2 A 

time (yrs) 



OH 




HDO 



0.01 



0.00 



gram 30K; G„ 10^ 




1 2 

time (yrs) 



H,0 



HDO 




timic (yrs) 



Fig. 2. Same as Fig. [T] with a UV radiation field of Go=10^. UV photons dissociate the chemical species present on the surface, 
which explains the low coverage on the grain surface (left panels). On the other hand, the amount of species released into the gas 
phase can be enhanced. 
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0.01 



0.001 



u 



10 



H30 



20 40 60 80 

time (years) 



100 




40 60 

time (years) 



100 



Fig. 3. Chemical species present on dust grains and released into the gas phase in diffuse clouds. Left panel: chemical species present 
on a grain surface in monolayer (1 monolayer is a total coverage of the surface) as function of time. Right panel: fraction of newly 
formed oxygen bearing molecules f(0) that are released into the gas phase under the form of other species. The error bars have been 
derived using 7 simulations and represent 95 % of confidence level. 



0.01 - 



>. 



0.001 



H3O 



1 20 30 

tim.e (years) 

Fig. 4. Same as Fig. |3]for translucent clouds. Note that no oxygen bearing species are released into the gas phase. The error bars 
have been derived using 4 simulations and represent 95 % of confidence level. 
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H,0 



u 



10 20 30 40 

time (years) 



50 




10 20 30 40 

time (years) 



50 



Fig. 5. Same as Fig. [3] for dense clouds. The error bars have been derived using 5 simulations and represent 95 % of confidence 
level. 






ILO 



/|0 60 

time (years) 



100 




0,001 



40 60 

tim.e (years) 



100 



Fig. 6. Same as Fig. |3]for PDRs. The error bars have been derived using 7 simulations and represent 95 % of confidence level. 
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Appendix A: Photo-desorption enhanced by photo-dissociation: recombination probability 

First, we assume that a photo-dissociation occurs, leading to one product on the (x = Q,y - 0) site and the other next to it. The 
probability Px,y(i), at a given position (x, y) on the grid, to find the product of the photo-dissociation after it has made / steps can be 
found through previous probabilities at the step / - 1 : 



Px.y = - [P,-l,y(i - 1) + Px+l,y(' " 1) + P^,y-l(' " 1) + ^^..v+l(' " 1)] 



(A.l) 



After further investigation, one can show that for the special position (x - 1,^ = 0), just next to the initial location, the probability 
can be written by separating even and odd cases: 



Pifii^i - 1) = 



(/ + If 



(A.2) 
(A.3) 



of 



Pifli^i) = 
More generally, the probablity Pij(n) that the product arrives after exactly n steps on a specified site situated at (x, y) is given by: 

n-l 

Px,y{n) (A.4) 



P{y{n) 






f f 

\fn-2k,keW, then one can conclude that Px,y{n) - 0, whereas \f n - 2k - \,k eW , Px,y(n) exists and can be simplified to: 



P{y(2k-1)^ 



k-l 



Y]{l-Px,y(2j-1)) 



Px,y(2k - 1). 



Hence for the case where x = 1 and y -0: 



P{/2k-l) 



n(i-4-'/) 



(k + iy 



In this last equation the product term can be seen as a function of integers (p(k) having for first terms: 



(A.5) 



(A.6) 



-'.•2 



1=2 
^(1) = 1 

>pi2) = 0.75 

.^(3) ^ 0.6445 

<p(4) ^ 0.6042 

lim (p(k) ^ 0.581739 



(A.7) 

(A.8) 

(A.9) 

(A. 10) 

(A.ll) 

(A. 12) 



Here we see how fast this function converges to an unique value. At the end, the probability Pifl that the initial product comes back 
to its site of production is then: 



Pi,o^Yu^(i)PiQ(i), 



(A. 13) 



where S (i) is the cumulative distribution of the number of steps that the product will do. According to the previous results that we 
obtained, this probability can be re-written as: 



Pi,o^J]S(2k-lMk) 



(k+lf 



k=l 



4/:+! " 

First, if we assume that the product can make as many steps as it wants (Vfc, S (2k - 1) = 1), Pio can be approximated as: 

'^ (k+lf 



k=5 



4*+! 



(A. 14) 



(A. 15) 



where cp* is the limit value when k -^ +00. To have a better idea of how fast this serie converges, note that the sum in the last term 
of the previous equation has an algebraically defined value which is very small compared to the other terms. Indeed: 



n 1 \2 

^•yillil ^0.581739.. 

■^ Zj 4k+i 27648 



365 
X r^— TT ^ 0.00768. 



(A. 16) 



k=5 
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Here we see that the direct recombination effect is enhanced mainly by the first term of the sum, i.e. when k < 5, hence a number of 
steps smaller than 9. This shows that one can assume that a product has left its initial position after 10 steps if it has not come back 
at least one time to its initial position (or neighbourhood here). Finally, we can estimate Pio to: 

Pi 0^0.4474, (A. 17) 

giving a new intrinsic probability of photo-desorption Pspd enhanced compared to the initial value. Indeed, if we take into account 
that a new molecule that is not photo-desorbed can be splitted by photo-dissociation, then the product can encounter for a second time 
the reactant, giving him a second chance to photo-desorb, and so on. Thus, an enhanced value of the spontaneous photo-desorption 
P* ^ can be expressed as: 

^spd = -Pspd + (1 - P.pd)PlflP.pd + (1 - P.pdPl,QfP.pA + ... (A.18) 

/■=() 

The last term being a geometric series that converges because |(1 - Pspd)^i,ol < 1, and one can simplify to: 

p; ^ . (A.20) 

^P^ l-(l-Pspd)Pl,0 

In the case of water formation with the reaction Oh- H —> OH and OH -i- H — » H2O, the fraction of water that desorbs upon formation 
is 5%. In the presence of a UV field, this value becomes 9.4%. 
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